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Abstract. We present the results of the broad-band photome- 
try of the nearby middle-aged radio pulsar PSR B0656+14 and 
its neighborhood obtained with the 6-meter telescope of the 
SAO RAS and with the Hubble Space Telescope. The broad- 
band spectral flux F v of the pulsar decreases with increasing 
frequency in the near-IR range and increases with frequency 
in the near-UV range. The increase towards UV can be natu- 
rally interpreted as the Rayleigh- Jeans tail of the soft thermal 
component of the X-ray spectrum emitted from the surface of 
the cooling neutron star. Continuation of the power-law com- 
ponent, which dominates in the high-energy tail of the X-ray 
spectrum, to the IR-optical-UV frequencies is consistent with 
the observed fluxes. This suggests that the non-thermal pulsar 
radiation may be of the same origin in a broad frequency range 
from IR to hard X-rays. We also studied 4 objects detected in 
the pulsar's 5" neighborhood. 
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1. Introduction 

The radio pulsar PSR B0656+14 was discovered by 

Its basic parameters are pre- 



Manchester et al (1978) 
sented in Table |l[ This pulsar has been detected in X-rays 
at E ~ 0.1-10 keV, with the space observatories Ein 

ROSAT 



Anderson et al. 



stein (Cordova et al. 1989) 



Foster et al. 199 



1993|), 
and 



EUVE 
ASCA ( jreiveldinger et al 



( Finleyetal. 1992 ; 
(Finleyetal. 1994; 



1996 



Zavlin et al. 200Cp . Together with Geminga, PSR B1055-52, 
and, perhaps, the Vela pulsar, it belongs to a group of 
middle-aged pulsars with thermal-like soft X-ray spectra 
(Becker & Triimper 1997, and references therein). The X-ray 



radiation of PSR B0656+14 can be described as a combination 
of radiation from the entire surface of a cooling neutron star 
(NS) and radiation from hotter regions, presumably polar caps 
heated by relativistic particles produced in the pulsar magneto- 
sphere. An excess over the hot thermal component at energies 
above « 2 keV was inte rpreted as nonthermal rad iation from 
pulsar's magnetosphere ( jreiveldinger et al. 1996 ). The pulsar 
has also been margi nally detected in gamma-r ays (E > 50 
MeV), at a 3<r level ( feamanamurthy et al. 1996 ). 

The optical count erpart of PSR B0656+14 was detected by 
|Caraveo et al. (1994)| with the ESO 3.6-meter and NTT tele- 
scopes in the V band (V» 25). It was observed in the long- 
pass filter F130LP (AA = 2310-4530 A) with the Faint Object 
Camera (FOC) aboard the Hu bble Space Telescope (HST) by 
Pavlov et al. (1996; hereafter P96), who showed that the bulk 



of radiation in this band is of a nonthermal origin. Pavlov et 
al. (1997; hereafter P97 ) measured the pulsar's fluxes in nar- 
rower FOC bands, F430W, F342W and F195W, and estimated 
thermal and nonthermal contributions to the optical-UV flux. 



Mignani et al. (1997) measured the flux in the F555W filter of 



the HST Wide Field Planetary Camera (WFPC2). Coherent op- 
tical pulsations with the radio pulsar period have been detected 
in the B-band by |Shearer et al. (1997)[ Kurt et al. (1998; here- 
after K98) presented first results of the ground-based BVRJ 
photometry of the pulsar with the 6-meter Big Alt-azimuth 
Telescope (BTA) of the Special Astrophysical Observatory of 
the Russian Academy of Sciences (SAO RAS). 

The optical-UV data obtained in the above-described ob- 
servations have shown that the pulsar's optical flux exceeds 
the Rayleigh-Jeans extrapolation of the thermal spectrum 



seen in soft X-rays. P97 fit the observed optical-UV fluxes 



Send offprint requests to: A.B. Koptsevich 



with a two-component model, a sum of the power-law and 
Rayleigh-Jeans spectra. The slope of the nonthermal compo- 
nent a p hot = — 2.4 ± 0.7, although poorly constrained in these 
fits, looks steeper than that in X-rays. The thermal compo- 
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Table 1. Parameters of PSR B0656+14 (Taylor et al. 1993). 



P 

ms 



Observed 



Derived 



P 
10 -i4 



DM 



pc 



/ 
deg 



b 
deg 



mas yr~ 



mas yr~ 



T 
Myr 



B 

G 



E 

erg s~ 



d b 
kpc 



384.87 5.50 



14.0 



201.1 8.3 



"Pulsar proper motion ( Vlignani et al. 2000) 

fir. .■ ... , ^ — - — : : — ~- r. 



43 ±2 



-2±3 



0.11 4.7 x 10 1 



3.8 x 10 d 



0.2-0.8 



a uuui uiuui/i iinjuwii ^ l *xiciiuiii vt ui. i-vy\jv_/ /. 

Estimates obtained from interpretations of X-ray data give d — 0.2-0.5 kpc; the dispersion-measure distance is ~ 0.76 kpc. 



nent is characterized by the Rayleigh- Jeans parameter, G = 
Ts(Rio / d^oo) 2 = 1-6, where Tq is the brightness tempera- 
ture in 10 6 K, R = 10 i?io km the apparent NS radius, and 
d = 500g?5oo pc the distance. To tighten the constraints on 
the parameters of the thermal and nonthermal components, the 
fluxes should be measured in a broader wavelength range. 

The pulsar field has been recently observed with the HST 
Near Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) in three broad bands, Fl 10W, F160W, and F187W (Har- 
low et al. 1998). The pulsar was detected in all the three bands, 
but the near-IR fluxes appeared to be significantl y low er than 
the extrapolation of the BVRI fluxes measured by £98. 



to as BVRI hereafter). Other technical details can be found in 



K98| also detected four unresolved, faint (V > 25) optical 
objects within 5" from the pulsar counterpart. Because of low 
spatial resolution of the ground-based observations, it was im- 
possible to determine whether they are just background objects 
(point-like or extended) or some structures associated with the 
pulsar. 

In this paper we report the results of new optical observa- 
tions of PSR B0656+14 and the objects in its nearest neighbor- 
hood taken with BTA. We incorporate the available HST NIC- 
MOS, WFPC2, and FOC observations^ analyze them together 
with the new and old BTA data, and obtain new estimates of 
the pulsar magnitudes in the BVRI and near-IR bands. We also 
discuss the morphology, color spectra and possible nature of 
the objects in the pulsar vicinity. 

The observations and data reduction are described in 
Sect. |J Multicolor photometry of the pulsar and the objects 
in its vicinity is presented in Sect. g. The broad-band spectra 
and the nature of the objects are discussed in Sect. |]. 



2. Observations and data analysis 

2.1. BTA observations 

We include in our analysis the data obtained in two BTA obser- 
vations - of 1996, November 11/12 (described by K98) and 
1997, November 26/27. Both observations were carried out 
with a CCD photometer, with a pixel size of 0"274 x 0'.'274. 
The response curves of the filters used are close to those of the 
B, V, R c , and I c filters of the Johnson-Cousins system (referred 



Observations with the NASA/ESA Hubble Space Telescope were 
obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Incor- 
porated, under NASA contract NAS5-26555. 



K98 



The 1997 observations were performed with the R and I 
filters only. Because of poor seeing and high sky background, 
we had to exclude from consideration seven of ten images taken 
with the I filter and three of five images obtained with the R 
filter. Observational conditions for the accepted exposures are 
presented in Table g. 

We reduced the data using the ESO MIDAS and STScI 
I RAF software. We summed the new R and I images with 
the corresponding images taken on November 1996. This de- 
creases the statistical errors of flux measurements in compari- 
son with the previous measurements based on 1996 data alone. 

The instrumental stellar magnitudes m mstr for each filter 
(B,V,R,I) are 

£ \ . k 

t, 



_instr 



-2.5 log 



j cxp 



5m 



cosZ 



(1) 



where / is the flux in counts (DN) for a given aperture, 5m 
is the correction for finite aperture derived from the point 
spread function (PSF) of bright stars, k is the extinction fac- 
tor, for which we used average values of the BTA observatory: 
k = 0.34, 0.21, 0.15, and 0.10 for the B, V, R, and I filters, re- 



spectively (Neizvestny 1983). The signal-to-noise ratios S/N 
and the magnitude uncertainties Am lnstr were calculated as 



N 



Ai 



= .f 



f 

G 


+ A 
2.5 


N 


[i + 


A 

^sky 



-1/2 



(2) 



(3) 



In 10 S ' 

where 07 is the standard deviation of flux in counts, A is the 
number of pixels in the source aperture, G is the gain, and A s k y 
is the number of pixels used for background measurement. The 
instrumental magnitudes b,v,r,i were then transformed to the 
Johnson-Cousins magnitudes B,V,R,I using the equations 



73-6 = 26.16(3) 

V-v = 26.28(2) 

R - r = 26.59(2) 

I -i = 25.67(2) 



• 0.08(3) -(6 - 
0.11(3) -(6- 
0.02(1)- (v 
0.06(2) • (r- 



(4) 



where the transformation coefficients were determined with the 



aid of secondary photometric standards described by K98. The 



absolute fluxes F (in erg cm 2 s l Hz 1 ) were calculated us- 
ing the following equation 

log J F = -0.4 (m + m°), (5) 
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Fig. 1. Le/f: Fragments of the BVRI 
images of the PSR B0656+14 field 
taken with the BTA during two nights, 
in 1996 and 1997. The images are 
smoothed over the area of 3 x 3 pix- 
els. The fragment sizes are 14"xl4". 
PSR B0656+14 (PSR) and the objects 
in its vicinity (o2-o5) are marked in 
the R image. Right: The correspond- 
ing contour maps. Crosses ("x") mark 
the PSR and o2-o5 objects centered at 
their positions in the HST WFPC2 im- 
age. Pluses ("+") mark the positions 
obtained by the center/gauss fit of the 
point-like objects on the BTA images 
(see Sect. 3.3 for details). The origin 
of the reference frame is placed at the 
WFPC2 pulsar's "x "-position; the axes 
are scaled in arcseconds. 
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Table 2. Log of the PSR B0656+14 observations with the BTA. 



Band 


Date 


t b 

''exp 
S 


'obs 

UT 


Z d 
deg 


Seeing 
arc sec 


Sky 
DN s _1 pix" 1 


Sky 
mag arc sec -2 


B 


11.11.96° 


2400 


- 


34.30 


1.7 


2.47 


21.96 


V 


11.11.96" 


2400 


- 


33.76 


1.5 


6.43 


21.20 


R 


11.11.96" 


3000 


- 


33.82 


1.5 


18.17 


20.44 




26.11.97 


500 
500 


00:49 
01:11 


31.82 
33.92 


1.6 

1.6 


20.91 
21.70 


20.30 
20.25 




Sum e 


4000 


- 


33.59 


1.6 


18.96 


20.40 


I 


11.11.96" 


3300 


- 


33.65 


1.4 


29.93 


19.04 




26.11.97 


600 
600 
600 


22:37 
22:49 
00:37 


33.56 
32.42 
30.90 


1.5 
1.5 
1.6 


29.32 
28.68 
30.46 


19.07 
19.09 
19.03 




Sum e 


5100 


- 


33.14 


1.5 


29.77 


19.05 



Summed data sets 



of 1996 (K98) 



Exposure duration 
Starting UT of the exposure 



Zenith distance 

Summed R and I exposures 

of 1996 and 1997 



with the zero-points provided by Fukugita et al. (1995)| 



m° B = 48.490, m° v = 48.613, 
m% = 48.800, m°j = 49.058. 



(6) 



2.2. HSTdata 

We also stu died the PSR B0656+ 14 field using the HST data 
obtaine d by |Mignani et al. (1997)| in the WFPC2 F555W band 
and by |Harlow et al. (1998)| in the NICMOS F110W, F160W 
and F187W bands. The standard reduction scheme described 
in the HST Data Handbook was applied to the WFPC2 im- 
ages processed by pipeline. The data were processed with the 
STScI I RAF software, summed with the I RAF task crrej, and 
a ~ 116" x 116" mosaic image containing three WFC and 
one PC field, corrected for the HST geometric distortion, was 
constructed. We applied the photometric formulae similar to 
Eqs. (|l|) and (g), with allowance for the WFPC2 cooldown and 
decontamination corrections. The correction for the finite aper- 
ture was performed with the aid of the PSF image gb71 61 3du 
from the WFPC2 PSF Library. 

For the primary reduction of the NICMOS data, we used 
the standard procedure calnicaQ For rejection of anomalously 
"cold" pixels, we applied the I RAF task cosmicrays to in- 
verted images. To remove the pedestal effectjj and to subtract 
the background, the I RAF task pedskyf] was implemented, 
with the same reference flatfield. It significantly suppressed 
the pedestal level, but increased the background level near the 



2 We made use of the calibration files i7M2297n, h421621an, 
h4214599n, i7H653gn, h2413242n, h241323pn, h241323qn, 
icbl 5080n, icbl 507pn, and icbl 507qn, and the dark calibration file 
created by Brian Monroe (NICMOS Helpdesk at STScI). 



See ittp://www.stsci.edu/instruments/nicmos/anom_pedestal.html 



See STAN (Space Telescope Analysis Newsletter) No. 18. 
www.stsci.edu/instruments/nicmos/NICMOS_stans/stan_018.asci 



edges of the images, which is not critical for our purposes since 
the objects of interest are near the image centers. To combine 
the individual dither images and get rid of cosmic rays, we 
made use of the package ditherll by ^ruchteret al. (1998) and 
followed the recipe described in the example supplied with the 
package releaseg. The shifts of single images were calculated 
making use of three bright stars in the field. The PSF images 
were generated by the Tiny Tim codeQ Aperture sizes were cho- 
sen to maximize the signal-to-noise ratio. 

To estimate the flux errors for the sources observed with 
the HST, we neglected possible systematic uncertainties which 
are smaller than statistical errors for faint objects, including the 
pulsar. Conversion of the fluxes to the Vega system magnitudes 
for the HST bands were performed using Eq. (|j) and the fol- 
lowing zero-points, calculated by the pipeline: 



^° -48.610, m? in = 49.305, 



'555 



110 



11 - 49.885, m? s7 = 50.148. 



(7) 



'160 



'187 



2.3. Astrometry 

For astrometrical referencing of the BTA images, we used the 
mosaic WFPC2/F555W image which covers a fraction of the 
BTA field. The coordinates of five relatively bright stars in this 
overlapping region were determined in the WFPC2 and BTA 
images. Then, a six-parameter coordinate transformation was 
constructed. The uncertainties of the transformation are equal 
to w 0'/07 for the R,I images and 0"15 (~ 0.5 pixel) for the 
B,V images. This allows us to relate the BTA and HST images 
with accuracy better than the BTA pixel size. 



ittp://www.stsci.edu/^fruchter/dither/di theiil.ps 



ittp://www.stsci.edu/software/tinytirn} 
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Fig. 2. WFPC2 image of the pulsar and surrounding field. The fragment is 1 l"6 x 1 1"6 in size and s moothed with FWHM = 
f/.' 1 8. The arrow in the left panel shows the direction of the pulsar proper motion ( Mignani et al. 2000|), its length corresponds to 
the distance traveled by the pulsar in 20 years. The rectangle and the asterisk are discussed in Sect. 3.4 



3. Multicolor photometry of PSR B0656+14 and objects in 
its neighborhood 

3.1. Morphology of PSR B0656+ 14 field in different bands 

14"xl4" fragments of the smoothed BTA BVRI imagesF 
containing the pulsar are presented in Fig. Ill In Figs. E 
and g we show a smoothed ll'.'6xll'.'6 fragment of the 
HST/WFPC2 image and non-smoothed 10" x 10" fragments of 
the HST/NICMOS images, respectively. All the images are ori- 
ented as indicated in Fig. |[ The contour plots of the images in 
Figs. |l|, H and [| are smoothed over the areas of 3 x 3, 7 x 7, and 
5x5 pixels for the BTA, WFPC2, and NICMOS images, re- 
spectively. The isophotes correspond to the following levels l n 
(in counts) above the background: 



I, 



S + nka, 



(8) 



where S is the mean sky value near the pulsar, n = 0, 1, . . . 
N; a is the sky standard deviation related to one pixel; k is the 
scaling factor. We used k = 1.25, 2 and 3,N = 10, 10 and 18, for 
the BTA, WFPC2 and NICMOS cont our p lots, respectively. All 
the images confirm the detection by K98 of five faint objects, 
ol-o5 (marked in Figs, fi] and 0), within ~ 5" around the pulsar 
position. 

Because of worse seeing during the BTA observations, we 
cannot determine from the BTA images whether the faint ob- 
jects are point-like or extended, except for o2 which is clearly 
extended and looks as associated with the pulsar. With the aid 

7 These and other images are available in FIT S format at 



of the WFPC2 and NICMOS images, we infer that the pulsar 
counterpart ol (a 20 oo = 6 ,l 59 m 48?19, <5 2000 = 14°14'21"3, as 
determined from the NICMOS imagesrh and the stellar object 
o5 (a 20 oo = 6 h 59 m 47!96, S2000 = 14°14'21'/9; ~ 3" west of 
the pulsar) have point-like source profiles, whereas o2, o3, and 
o4 are extended. 

The extended sources are most clearly resolved in the NIC- 
MOS/F160W band. The object o2 shows a compact bright nu- 
cleus at a 20 oo = 6 h 59 m 48! 16, <5 20 oo = 14°14'23"2 (1"85 north 
of the pulsar). The nucleus is surrounded by a faint elongated 
nebula which occupies an area of ~ 3" x 0"7. The minimum 
distance between the nebula's southern edge and the pulsar, 
~ l'.'l in the F160W band, is smaller than the BTA seeing. 
This explains naturally the apparent association of the pulsar 
and o2 in the R and I bands. In the F555W band, whose re- 
sponse curve is very close to that of the Johnson's V filter, o2 
looks like a faint elongated clump without a pronounced nu- 
cleus. Also, o2 is barely seen as a very faint extended source in 



the FOC/F130LP image obtained by|P96. 

The extended object o3 (a 2000 = 6 h 59 m 48!23, <5 200 o = 
14°14'25"1) is ~ 4" north of the pulsar. It has an irregular 
structure with a bright knot at its eastern edge. Being at only 
~ 1" from the northern edge of o2, the object o3 appears to be 
associated with o2 in the F160W image. The structure of the 
object o4 (a 2000 = 6 /l 59 m 47!84, S 200 q = 14°14'22"0; - 5" 
west of the pulsar) varies strongly from band to band. 



http://www.ioffe.rssi.ru/astro/NSG/obs/0656-phot.html 



8 This and other absolute positions in the HST images are given as 
directly measured with the I RAF imcntr task, without correcting for 
possible HST pointing errors. 
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Fig. 3. The same as in Figs. |J g but for the HST/NICMOS bands. The fragments represent the region 10" x 10" in size of the 
non-smoothed images. 
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Table 3. Flux measurements of the pulsar's optical counterpart. 



Band 


S/N 


r a 


Source 6 


k/cos{Z) c 


"fin 


Mag e 


Flux 






pix 


DNs" 1 


mag 






/Jy 


B+ 


10 


3 


0.839(80) 


0.41 


2.44 


24.98(13) 


0.409(36) 


V+ 


9.0 


3 


1.26(14) 


0.25 


2.07 


24.97(12) 


0.369(39) 


R+ 


16 


3 


2.77(17) 


0.18 


2.15 


24.48(7) 


0.488(32) 


1+ 


10 


3 


2.26(22) 


0.12 


2.05 


23.87(11) 


0.675(65) 


B x 


8.1 


3 


0.661(81) 


0.41 


2.44 


25.25(13) 


0.319(36) 


V x 


9.2 


3 


1.25(14) 


0.25 


2.07 


24.97(12) 


0.369(39) 


R x 


11 


3 


2.24(19) 


0.18 


2.15 


24.71(9) 


0.394(31) 


I x 


9.5 


3 


2.02(21) 


0.12 


2.05 


23.99(11) 


0.604(58) 


F555W 


23 


4 


0.0853(38) 


- 


1.12 


24.90(7) 


0.393(24) 


F110W 


11 


2 


0.119(11) 


- 


1.55 


24.38(10) 


0.336(30) 


F160W 


14 


3 


0.185(13) 


- 


1.50 


23.22(8) 


0.575(41) 


F187W 


8.0 


3 


0.127(16) 


- 


1.65 


22.62(13) 


0.779(97) 



Aperture centered at the position determined by the center/gauss fit in the BTA images 

Aperture centered at the positions determined from the WFPC2 image, (see '+' and ' x ' marks in Fig. hh 

Aperture radius in pixels 

Source count rate within the aperture 

Correction for atmospheric extinction 

Factor to correct for the finite aperture 

Johnson-Cousins and Vega magnitudes, for the BTA and HST images, respectively 



3.2. Consistency of the BTA and HST data 

To check consistency of our data reduction and photometry of 
the objects observed with the different telescopes and detectors, 
we measured the magnitudes of a relat ively bright star (marked 
N7 in Fig. g, following the notation by frC98 ) seen in all the BTA 
and HST images. The raw broad-band spectrum of this star 
(Fig. |] and Table Q) looks rather smooth and shows no jumps 
between the BTA and HST bands. Moreover, the fluxes in the 
almost equivalent V and F555W bands are in good agreement 
with each other. This indicates that there are no major inconsis- 
tencies between the reduced BTA and HST data. The colors and 



magnitudes of N7 correspond to an M5V star (Bessel 1990 1 at 
a distance of 440 pc. 

3.3. Flux measurements 

The objects of interest are faint and strongly blurred in the 
BTA images, especially in the R and I filters. This leads to 
additional uncertainties of their center positions and fluxes. To 
evaluate the uncertainties, we determined the positions (and the 
corresponding fluxes) of the unresolved objects using two ap- 
proaches. First, we applied the center/gauss task of the MI- 
DAS package to our BTA images. These positions are marked 
by the '+' signs at the contour plots in Fig. fil (The center of 
o5 cannot be determined by this method in the B and V im- 
ages because o5 is hidden in the background.) Second, we mea- 
sured the positions of ol (PSR), o5 and o3 in the WFPC2 im- 
age and transformed them to the BTA positions as described 
in the Sect. 2.3. These positions are marked by the 'x' signs 



(o4 cannot be centered in this way because it is resolved in the 
WFPC2 image.) We see from Fig. [y that the '+' positions of 
the point sources (ol and o5) are systematically closer, by 1-2 



BTA pixels, to the nearest extended objects (o2 and o4, respec- 
tively) than the ' x' positions. This means that the fluxes of the 
faint, point-like sources may be contaminated by the nearby ex- 
tended objects, and the Gaussian fits of their profiles may not 
be reliable. 

We measured the stellar magnitudes and fluxes of point- 
like objects for a range of aperture radii r. For each of these 
objects in the BTA images, we chose r = 3 pixels as an optimal 
value and measured the background over the annulus with inner 
and outer radii r- m = 8 and r ou t = 11 pixels, centered at the 
source. To measure the fluxes in the HST images, we used the 
task phot from the I RAF package apphot. For the NICMOS 
images (pixel size is / .'075), we chose r — 3 pixels (except 
for the pulsar in Fl 10W - see Table pb, and r- m — 6, r out = 9 
pixels (ri n = 12, r out = 15 pixels for N7). We used r = 4, 
r in = 15, r out = 20 pixels for the WFPC2 image (pixel size is 
(X'045). The results are presented in Tables andH. 

We see from Table [3] that the B, R and I fluxes of the pul- 
sar, measured with the apertures centered on the ' x ' positions 
are systematically lower, by ~ 10-30%, than those obtained 
with the '+' positions. In the V band both fluxes coincide with 
each other and are consistent with the F555W flux, within the 
errors. Similar differences between the 'x' and '+' fluxes are 
found for o3 if it is considered as a point-like source. Further- 
more, the ' x ' fluxes yield a smooth broad-band spectrum of o5 
(see Table and Fig. Eh while its '+' flux in the I band stands 
out of the smooth spectrum by ~ 30%, an excess hardly plau- 
sible for an ordinary object. Based on these measurements, we 
believe that the ' x ' fluxes of the point-like sources are more 
realistic than the '+' fluxes, and even the ' x ' fluxes are likely 
overestimated because of contamination from nearby extended 
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Table 4. Fluxes (upper values, in /iJy) and Johnson-Cousins (BVRI) or Vega magnitudes (lower values) of the objects in the 
vicinity of PSR B0656+ 14. 



Object 


B 


V 


F555W 


R 


I 


F110W 


F160W 


F187W 


o2 


0.122(34) 
26 m 29(30) 


0.187(38) 
25™71(22) 


0.36(12) 
25 m 00(36) 


0.211(32) 
25 m 39(16) 


0.526(59) 
24 m 14(12) 


1.14(14) 
23 m 06(13) 


3.58(15) 
21 m 23(4) 


4.93(27) 
20 m 62(6) 


o3 


0.205(36) 
25 m 73(19) 


0.282(39) 
25™26(15) 


0.35(11) 

25 m 03(34) 


0.353(34) 
24^83(10) 


0.512(62) 
24™ 17(13) 


0.77(16) 
23™49(23) 


1.66(15) 
22™ 07(10) 


2.40(27) 
21™ 04(1 2) 


o4 


0.217(38) 
25 m 67(19) 


0.136(34) 
26 m 05(27) 


0.42(12) 
24 m 83(31) 


0.398(32) 
24™ 70(9) 


1.099(69) 

23™ 34(7) 


1.00(13) 
23™ 19(14) 


1.64(12) 
22 m 08(8) 


1.44(25) 
21 m 95(19) 


o5 


- 


- 


0.268(20) 
25 m 32(8) 


0.520(41) 
24 m 41(9) 


0.774(75) 
23 m 72(ll) 


0.991(33) 
23 m 21(4) 


1.488(37) 
22 m 18(3) 


1.374(78) 
22 m 01(6) 


N7 


6.37(34) 
22 m 00(6) 


26.7(12) 
2(^32(5) 


23.364(40) 
20 m 469(2) 


64.9(12) 

24 m 17(2) 


192.8(52) 
17 m 73(3) 


313.99(25) 
16 m 954(1) 


439.96(25) 
16 m 008(l) 


386.65(30) 
15 m 884(1) 



We also estimated integral fluxes of the extended sources 
in different bands. For consistency, the fluxes were measured 
with the same circular aperture for all the objects and bands; its 
diameter, l'/64, was chosen in accordance with the BTA see- 
ing conditions in the I band. To measure the fluxes of o2, we 
centered the aperture at the center of its nucleus in the F160W 
image (this position was transformed into the BTA and WFPC2 
images using the pulsar and o5 as reference objects). Although 
this aperture does not include the whole o2 nebula, the contri- 
bution of the outer parts of the nebula to the integral flux is only 
a few percents, smaller than the errors of our measurements. 
To estimate the BTA BVRI fluxes, we corrected the measured 
fluxes for the finite aperture as if o2 were a point-like object. 
Such corrections are negligible for the HST images. The same 
approach was used for o3 and o4. For these objects, the contri- 
butions into the measured integral fluxes from the nearby ob- 
jects o2 and o5 are not significant. 

The magnitudes and the raw broad-band spectra of the ex- 
tended objects are presented in Table |] and Fig. |J For o2 and 
o3, the spectra are similar and can be crudely approximated by 
a power law, F v oc v~ a , with a « 2.4 and « 1.6, respectively. 
The spectrum of o4 appears to be more complicated. However, 
the estimated fluxes of this object are less definite because of 
its clumpy morphology and rather low S/N (~ 1.5-2), partic- 
ularly in the B, V and F555W bands. 

3.4. Contamination of the pulsar flux by the object o2 in the R 
and I bands 

We see from Tables and H that the ' x ' fluxes of the pulsar 
in the R and I bands are comparable with the corresponding 
fluxes of o2 whereas the '+' fluxes are higher in both filters. The 
reason is that the '+' positions are shifted towards o2 by ~ 0'.'4, 
so that the '+' fluxes are stronger contaminated. Therefore, we 
consider the ' x ' fluxes more reliable. 

If o2 had a point-like profile, its contribution in the pulsar's 
aperture would be negligible, about 3% in the I band and even 
less than that in the R band. If o2 had the same surface bright- 



ness at the pulsar position as the brightness measured over the 
rectangular aperture of 4.6 arcsec 2 (see Fig. 0), its contribution 
to the flux measured in the pulsar's aperture would be 84±7% 
and 96±6% in the R and I bands, respectively. Since o2 is not 
a point-like object, and it is spatially separated from the pul- 
sar, the obtained estimates can be considered only as lower and 
upper limits on the contamination. 

An intermediate, more realistic estimate can be obtained 
based on the fact that the brightness distribution of o2 in the 
F160W image is almost symmetrical with respect to its major 
and minor axes. We assumed the same source symmetry in the 
I and R bands, and measured the I and R fluxes in the aperture 
with the radius of (X'82 (3 pixels) centered at the point marked 
by the '*' sign in Fig. g, correcting the fluxes for the finite PSF 
as if there were a point source. No excess over the background 
level was detected at this position in the R image. In the I image 
we detected F* = 0.16 ± 0.06 /iJy, i.e., 27±10% of the pulsar 
flux F* , which can be considered as an estimate for the con- 
tamination. This value is comparable with the 3a uncertainty 
of the flux measured in the pulsar's aperture as well as with the 
difference between the Fj~ and Ff fluxes (see Table ph. No- 
tice, however, that some contribution to the flux at the '*' po- 
sition can also be provided by o3, so that the real contamination 
from o2 may be somewhat lower than estimated above. More- 
over, the assumption on symmetry of o2 in the I band should be 
taken with caution because, according to Fig. [lj the brightness 
distribution of o2 (and, perhaps, of o3) in this filter appears to 
be different from those in the other bands. 

The exposure in the B band was too short to estimate the 
contamination. Since Fg — Fg ~ 0.1 [i]y is comparable to 
Fj — Fp , we can expect that the pulsar's B flux can also be 
contaminated by o2. Although this object is not seen in the ad- 
jacent F430W filter and in the other HSTUV bands, its south- 
ern fragment is seen in the much wider F130LP filter whose 
band-pass strongly overlaps with that of the B filter. 

Given the large uncertainty of the contamination, it is very 
difficult to obtain accurate values of the B, R an d, particularly, 
I fluxes. Since the F + fluxes, measured by K98, are obviously 
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Broad-band spectra of the stellar objects o5 and N7. 



overestimated, we can choose either F x or F x — F* as most 
probable values. Since F^ <C F£, and Fg is difficult to mea- 
sure, we choose F£ and Fg as "best-guess values". As for the 
I flux, we present both the Ff and Ff — F* values in Table g, 
where we collected all optical data for the PSR B0656+14. 
Only observations with spatial resolution better than 1" can im- 
prove the accuracy of the flux measurements in the BRI bands. 



4. Discussion 

4.1. Spectrum of PSR B0656+ 14 

In Fig. ra we combined our optical data with the available mul- 



tiwavel ength data on PSR BQ6 56 +14 (fVIalofeev 1999J; |P96 ; 
?9l\ |Korpela & Bowyer 1998fc |Ramanamurthy et al. 1996; 



Zavlin et al. 200C). The spectrum generally fades towards 



higher frequencies, showing a hump of thermal radiation in 
soft X-rays. 

In X-rays, the EUVE, ROSAT and ASCA data fit well 



with a three-component spectral model (Zavlin et al. 2000). 
The model consists of a soft blackbody component describing 
thermal emission from the entire NS surface with a radius R s = 
13.5 ± 1.0 km and temperature T s = (0.84 ± 0.03) x 10 6 K, a 
hard blackbody component representing thermal emission from 
pulsar hot polar caps with T h = (1.65 ± 0.19) x 10 6 K and 
Rh = 0.9 ± 0.5 km, and a power-law non-thermal component, 
presumably formed in the NS magnetosphere, Je = CE a 
phot cm~ 2 s _1 keV -1 , where E is photon energy in keV, 
a = -1.45 ± 0.26 and C = (2.8 ± 0.1) x 10~ 5 . The radii 
R s and Rh are related to the pulsar distance d = 500 pc; 
column density of the interstellar matter towards the pulsar is 
uh = (1. 42±0. 06) x 10 20 cm~ 2 . The sum of the three compo- 
nents is displayed with the bold line in Fig. pi the unabsorbed 
spectrum is shown with the dotted line; the dashed line repre- 
sents an absorbed extension of this three-component model to 
the optical domain; the dot-dashed lines represent contributions 
of different components of the fit. The soft thermal component 
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from the entire surface of the cooling NS dominates in soft X- 
rays and in the EUVE band. The non-thermal component is 
responsible for the high energy tail in X-rays. It is also appar- 
ently detected in 7-rays, but its slope is highly uncertain and 
likely different from that in X-rays. 

The IR-optical-UV range is blown up in the inset of Fig. H. 
It is seen that the observed broad-band spectrum is consistent, 
within the uncertainties, with the continuation of the X-ray fit 
to the optical range. A typical uncertainty of the extrapolated 
flux is indicated by the dashed error-bar on the left side of the 
inset. The power-law component becomes dominating in the 
optical-IR range. This suggests that a single mechanism of non- 
thermal radiation may operate in a broad frequency range from 
IR to hard X-rays. The slope of th e non -ther mal component is 
less steep than it was deduced by K98 and P97 based on the 
previous optical-UV data. The overall optical spectrum has a 
shallow minimum at v « 7.4 x 10 14 Hz (A « 4000 A). Above 
this crossover frequency, the Rayleigh-Jeans part, F v oc v 2 , of 
the soft thermal component apparently dominates in the pulsar 
spectrum (P97). 

From the new BTA and NICMOS data we cannot exclude 
nonmonotonous behavior of the broad-band IR-optical spec- 
trum, which might indicate the presence of unresolved spec- 
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tral features. There are an apparent excess in B (K98; P97), 
an increase towards RI, and a dip in F110W, accompanied 
by a steep increase towards the IR frequencies. The BRI ex- 
cesses could be partially due to contamination of the pul- 
sar flux by the nearby extended object o2. In the R and B 
filters they could be caused by H Q and Hp emission from 
an unresolved pulsar-wind nebula (PWN) around the pulsar 



Table 5. Pulsar fluxes in the optical-UV bands. 



(e.g., Cordes 1996 1. For typical ISM temperatures and densi- 
ties, the pulsar's velocity, v p <; lOOc^oo km s _1 is super- 
sonic, so that the PWN would be a bow shock at a stand-off 
distance R w » (E/AircpuQ 1 / 2 = O.OOlZE^n-^v'l pc, 
where n — p/rtiH is the baryon number density of the ambi- 
ent medium, v-n.t = w p /(100 km s _1 ). For instance, such a 

= 2.5 x 10~ 3 photons cm~ 2 
S _1 , was observed at 7" (R w = 0.006 di 80 pc for £'( 0437 ) = 
1.2 x 10 34 erg s" 1 ) from the millisecond pulsar J0437-4715 



,J pJ 
shock, with an H a flux f H 



( |Bell et al. 1995) 1. Since f Ha ex Ev p Xd~ 2 , where X is the 
fraction of neutral hydrogen atoms, the expected H Q flux from 
a bow shock at ~ 2"-3" from PSR B0656+14 can be estimated 
by simple scaling: /^ 656) = 1.0 x i -3x( 0656 )/X( 0437 ) pho- 
tons cm" 



- 1 , fori; 



(0437) 



(0656) 



100 km s" 1 ,^ 0656 ) 



500 pc, d( 0437 ) = 180 pc. This corresponds to an R-band en- 
ergy flux i^ 0656) « 2I(° 656 )/I(° 437 ) pjy from the whole H a 
emitting region. The contribution of this flux in the observed R 
flux would be very substantial at X^ 0656 ) = X^ 0437 K How- 
ever, PSR B0656+14 is much younger than J0437-4715, and 
its hot NS creates an ionization zone of a radius R s ~ 0.3-0.5 
pc (~ 2'-3') around the pulsar (e.g., jlhvartsman 197d| ). Since 
X 1 - 0656 ) at r — R w (-C R s ) is expected to be orders of mag- 
nitude smaller than X^ 0437 ' , the contribution of the bow-shock 
H Q emission in Fr should be negligible. The estimated angular 
size of the ionization zone is much larger than our apertures, its 
emission is hidden in background and cannot contribute to the 
measured point source flux. 

A similar excess, albeit in the V band, is also present in 
the broad-band spectrum of the other well-studied middle- 
aged pulsar Geminga. Its possible explanations have been 



suggested by Bignami et al. (1996) and Martinet al. (1998) 



based on various combinations of the synchrotron or ion 
cyclotron emission/absorption in an upper atmosphere or 
low magnetosphere of t he NS. In Fig. \l we plot the 
optical-UV fluxes (from [Vlignani et al. 1998 ) of Geminga 
together with our estima tes of the NICMOS F110W and 
F160W fluxes detected by |Harlow et al. (1998)[ We also show 
the available dereddened color spectra of the Crab pulsar 
d Percival et al. 199$ [Eikenberry et al. 1997]), PSR B0540-6 9 



(Middleditchetal. 1987), the Vela pulsar ( Masuti et al. 1997) 



and PSR B0656+14 (see Table g). We neglected the extinc- 
tion for PSR B0656+14 and Geminga, which gives rather small 
corrections for these nearby pulsars. It is seen that the broad- 
band spectra of PSR B0656+14 and Geminga are qualitatively 
similar to each other. This suggests that the unresolved spec- 
tral features associated with the apparent optical excesses in 
their spectra could be of the same nature. On the other hand, 
the spectra of the middle-aged pulsars appear to be different 



Band 


M a 


Flux 


References 6 






^Jy 


data / measured by 


B 


25.251^ 


0.32iJ° 


tw 


V 


24.97(12) 


0.369(39) 


tw 


R 


74 77+25 


0.39±J° 


tw 


r 


23.99(11) 


0.604(58) 


tw 




24.33(23) 


0.441(83) 




F555W 


24.90(07) 


0.393(24) 


M97a / tw 


F110W 


24.38(10) 


0.336(30) 


H98 / tw 


F160W 


23.22(08) 


0.575(41) 


H98 / tw 


F187W 


22.63(13) 


0.779(97) 


H98 / tw 


F430W 


- 


0.264(29) 




P97 




F342W 


- 


0.307(27) 




P97 




F195W 


- 


0.349(43) 




P97 




F130LP 


- 


0.305(06) 




P96 





Johnson-Cousins (for the BTA images) or Vega system 
(for the HST images) magnitu de 



tw 



M97a - |Mignani et al. (1997)|, H98 -|Harlow et al. (1998) 
this work 



See explanations in Sect. 3.4 



from the flat and featureless spectra of the young Crab and 
B0540— 69 pulsars, which is evidence of evolution of optical 
emission with pulsar age. To confirm the difference and reveal 
its nature, spectroscopic observations of the middle-aged pul- 
sars are needed. 



4.2. Nearby objects 

Object o5. Faintness of the stellar object o5 suggests that it 
might be a main sequence star located near the edge of our 
galaxy. Taking into account the total interstellar galactic ex- 
tinction towards the pulsar field (I = 201?1, b = 8?26; A B = 
0.420, A v = 0.323, A R = 0.260, Aj_ = 0.189, A., = .088, 

Schlegel et al., 1998J ), we 



A H = 0.056, A K = 0.036 - see 

obtain the dereddened color indices of o5: Mp^f^ w — R 

0.85(9) and R - 1 = 0.62(12). This would be compatible with 



a K6V main sequence star (Bessel 1990). However, typical ab- 
solute magnitude of the K6V stars are My = 8 r . n 2-8 r . n 5. This 
yields a distance range of 20-23 kpc, incompatible with the 
Galaxy size. Therefore, it is more likely that o5 is a white dwarf 
with T - 4000 K and M v - 16 m , at a distance of - 600 pc. 
However, in the absence of spectral information one cannot rule 
out also that this is a quasar. 

Extended objects. There were unsuccessful attempts to 
search for a PWN around PSR B0656+14 on scales ^ 
10" in radio (Frail & Scharringhausen 1997) and X-rays 



( [Becker et al 1999[ ) 

{E/(47rcp )Y^ 



However, an expected PWN size, ~ 
where po is the pressure of the ambient 



medium, can be smaller than 10", and, in principle, the object 
o2 might be a PWN candidate. The morphology of o2 looks 
like that of a faint spiral galaxy seen edge-on. However, its 
size, ~ 3" x 077, is larger than the mean size, ~ 1" x 075, 



of galaxies of close magnitudes, 



l F160W 



22, in the South 
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Fig. 6. Broadband spectrum of PSR B0656+14. The inset blows up the optical range of the spectrum. 
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Hubble Deep Field NICMOS catalog| taken with the exposure 
4992 s close to the pulsar field exposure 5098 s in the F160W 
band. The broad-band spectrum of o2 is close to a power law 
with a photon index of about —3.4. If it were a nebula pow- 
ered by relativistic particles ejected by the pulsar, its observed 
flux could be provided by synchrotron radiation from a region 
with approximately the same parameters as those estimated by 
Hester et al. (1995) for the optical structures (wisps, knots) in 
the synchrotron nebula around the much more energetic Crab 
pulsar (e.g., from a volume V ~ 2 x 10 48 cm 3 with a mag- 
netic field B ~ 3 x 10~ 3 G and a steep power-law distri- 
bution of relativistic electrons with a slope p — 5.8, number 
density n e ~ 3 x 10' 5 cm' 3 , and minimum Lorentz factor 
7min ~ 10 4 ). However, the morphology of o2 looks very un- 
natural for a PWN, particularly with account of the pulsar's 
proper motion (see Fig. 2), the slope of the power-law spectrum 
is too steep compared with other known synchrotron sources, 



See nicmos_vl generic. cat, presented at 



and PSR B0656+14 can hardly supply relativistic particles at a 
rate needed to support the amount required. Thus, we conclude 
that o2 is most likely a faint background galaxy, as are the other 
extended objects o3 and o4. 

Deep subarcsecond optical, radio and X-ray observations of 
the PSR B0656+14 field are needed to search for possible un- 
resolved structures associated with activity of PSR B0656+14. 
In this respect, it would be useful to confirm and investigate a 
knot-like structure seen at a 3cr level above the background at 
~O704 from the pulsar position in the FOC/F342W imageQ- 
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